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ABSTRACT 

The generation of mean magnetic fields is studied for a simple non-helical flow where a net 
cross helicity of either sign can emerge. This flow, which is also known as the Archontis 
flow, is a generalization of the Arnold-Beltrami-Childress flow, but with the cosine terms 
omitted. The presence of cross helicity leads to a mean-field dynamo effect that is known as 
the Yoshizawa effect. Direct numerical simulations of such flows demonstrate the presence of 
magnetic fields on scales larger than the scale of the flow. Contrary to earlier expectations, 
the Yoshizawa effect is found to be proportional to the mean magnetic field and can therefore 
lead to its exponential instead of just linear amplification for magnetic Reynolds numbers that 
exceed a certain critical value. Unlike a effect dynamos, it is found that the Yoshizawa effect 
is not noticeably constrained by the presence of a conservation law. It is argued that this is 
due to the presence of a forcing term in the momentum equation which leads to a nonzero 
correlation with the magnetic field. Finally, the application to energy convergence in solar 
wind turbulence is discussed. 
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1 INTRODUCTION 

The dynamo effect in astrophysical objects is often associated with 
the occurrence of helicity in them. In magnetohydrodynamics there 
are several helicities that can be important. A particularly impor- 
tant one is the kinetic helicity, because its va lue is finite in r otating 
stratified bodies and can lead to an a effect iMoffatl|[T978l : |Parkeii 
Krause & Radleil ll980^. Another important helicity is the 
magnetic helicity. Unlike the kinetic helicity, the magnetic helic- 
ity is conserved by the quadratic interactions, so its value can only 
change through resistive effects or th rough magnetic helicity fluxes 
( [Brandenburg & Subramani^l2005l) . Such a conservation law is 
crucial to understanding the saturation behavior of a effect dy- 
namos. This is because the a effect tends to produce large-scale 
magnetic fields that are helical, but conservation of total magnetic 
helicity implies that there must be small-scale magnetic helicity of 
the opposite sign, so that the sum of small-scale and large-scale 
magnetic helicities is close to zero. This then leads to a resistively 
slow saturation phase in the nonlinear regime (Brandenburg 2001 ). 
Mathematically, the consequence of magnetic helicity conserva- 
tion can be described by the attenuation of the total a effect by 
the addition of a term proportional to the magnetic helici ty effect 
dField & BIackmarj|2002l : iBlackman & Brandenbur3l2002h . 

In a topological sense, magnetic heli city describes the link- 
age of magnetic flux tubes jMoffattlll969l) . while the kinetic he- 
licity characterizes the linkage of vorticity tubes. However, there 
is yet another helicity, the cross helicity, that describes the linkage 
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of magnetic flux tubes with vortex tubes. This quantity is impor- 
tant because it too is conserved by the quadratic interactions, i.e. 
it can change only by visco-resistive effects or by cross helicity 
fluxes. Moreover, the sma ll-scale cros s helicity can itself lead to 
large-scale dynamo action lYoshizawa : 1990). Such a mechanism is 
quite different from the a effect, because it corresponds to an inho- 
mogeneous term in the dynamo equations a nd could therefore play 
the ro le of a turbulent battery term. Indeed, [Brandenburg & UrpinI 
( Il998h showed that the battery term due to cross helicity can facil- 
itate large-scale dynamo action in young galaxies and hence could 
be responsible for the relatively strong magnetic fields observed in 
such galaxies at high redshifts. 

I n spite of several additional studies dYoshizawa & YokoH 

1 19931 : lYokoilll996l ; iBlackman & Choij|l997 '). large-scale dynamo 
action due to cross helicity has not received much attention be- 
cause this effect was never seen in simulations, nor was it found to 
be responsible for driving large-scale magnetic fields found therein. 
Such an effect would require that the small-scale magnetic field is 
systematically aligned with the flow, i.e. it is either mostly parallel 
or mostly anti-parallel to the flow. Such circumstances are known 
to prevail in the solar wind, but here the field comes presumably 
directly from the Sun and would therefore not be produced by a 
dynamo. 

In t he pr esent paper we consider the so-called 
Archontis! ilOOd) dynamo (see also iDorch & Archonti j |2004 



Cameron & Gallowavl[200^ which is driven by a forcing function 



that is based on the Amold-Beltrami-Childress (or ABC) flow, but 
wit h the cosine terms be i ng om itted. This flow was first proposed 
by iGallowav & Proctoil ( Il992h to study fast dynamo action by 
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calculating growth rates for the kinematic version of this flow. 
The ABC flow is helical and produces efficient dynamo action 
dOallowav & FrischI Il986h . However, the omission of cosine 
terms renders the flo w nonhelical, so that ther e is no a effect, 
but numerical studies toorch & Archontisll200i) have shown that 
such a dynamo produces magnetic fields that are either aligned 
or anti-aligned with the flow almost everywhere. This means that 
there is cross-helicity in the system, which can give rise to the 
Yoshizawa effect and produce large-scale dynamo action. Further- 
more, owing to the conservation property of cross helicity, such 
dynamos may be controlled by this effect and may also show slow 
saturation behavior. It is therefore of interest to investigate whether 
the formulation for the slow saturation of a effect dynamos carries 
over to the present case. 

We begin by explaining first the simulations, discuss the fea- 
tures of the kinematic growth phase of the dynamo, and then con- 
sider the slow saturation regime using a nonlinear dynamical feed- 
back formalism that is analogous to the dynamical quenching for- 
malism for the a effect. Next we argue that the kinematic growth 
in such a dynamo is indeed due to the cross helicity effect. We 
show that the estimated growth rate obtained from a simple model 
involving the induction and momentum equations along with the 
evolution equation for the small-scale cross helicity can be brought 
in good agreement with our simulation results. 



2 BASIC EQUATIONS 

We consider here a model that is similar to that o f lArchontislfcoOOh 
and iDorch & Archontij j2004) who assumed a compressible gas 
with an energy equation included. However, in their model the tem- 
perature was kept approximately constant by applying a heating 
and cooling term. Here we assume instead an isothermal equation 
of state, i.e. the pressure is given by p = pCg , where p is the density 
and Ca is the isothermal sound speed. The evolution equations for 
the density p, velocity U, and magnetic vector potential A are then 

2^ = -dV lnp + f+ i[JxB + V- i2puS)] , (2) 
Dt p 

^ = UxB + vV^A, (3) 
ot 

where D/Df — d/dt + U ■ ^ is the advective derivative, B — 
V X j4 is the magnetic field, J — V x B//io is the current density, 
fio is the vacuum permeability, 77 is the magnetic diffusivity, which 
is assumed constant, v is the kinematic viscosity. 



ko = Cs = po = Mo = 1, 



(6) 



S.i = |(C/.,i + C/,,0-|'5,,V.j7 
is the traceless rate of strain tensor, and 

F — Fo (sin koz, sin kox, sin koy) 



(4) 



(5) 



is the forcing function where Fo is an amplitude factor and ko is a 
wavenumber. 

For analytic considerations we consider the flow to be incom- 
pressible, i.e. V ■ [/ = and p — po = const. While this sim- 
plifies the treatment significantly, it should be remembered that the 
differences between compressible an d incompressible cases are no t 
critical if the Mach number is small dCameron & Gallowavlbooeh . 
In the present paper we consider cases where the Mach number is 
around 0.03 (see below). In order to simplify the notation we use 
units where 



although in several places we shall keep these units for clarity. 

The simulations have been performed using the PENCIL 
CodeB Triply periodic boundary conditions are employed for all 
variables over a cubic domain of size L x L x L. As initial con- 
dition we use zero velocity, constant density given by p = po, and 
a spatially random vector potential of sufficiently low amplitude so 
as to obtain a clear initial exponential growth phase over several or- 
ders of magnitude before nonlinear effects become important and 
lead to saturation of the magnetic field. 

Our simulations are characterized by the values of the mag- 
netic Reynolds and Prandtl numbers. 



R 



and Pm = 



(7) 



Up 

r/ko V 

respectively. Here, we have defined uo = {Fo/ko)^^^ as our 
reference velocity. Occasionally we also use the visco-resistive 
Reynolds number, 

R 



Ru — 



no 
pko 



1 



(8) 



where p — 1/ + rj. Throughout this paper we restrict ourselves to 
the case Pm = 1. The forcing amplitude is chosen such that the 
Mach number. Ma = uq/cs is small (about 0.03), so the flow stays 
close to incompressible. 

The flow is of course isotropic with respect to the three coordi- 
nate directions, so there is no preferred definition for the mean field 
in this case. Indeed, there are three equivalent definitions of two- 
dimensional averages (xy, yz, and xz averages). They all would 
lead to finite mean flows and mean magnetic fields. In the follow- 
ing we consider mean fields defined by averaging over the x and y 
directions, i.e. 



B{z,t) 



B dx dy. 



(9) 



Throughout this paper we focus on the case L = Lo, where we 
have defined Lo = 2n/ko. However, on one occasion we compare 
with the cases L — 2Lo and 4Lo, where the domain is big enough 
to allow for a field configuration that is four times bigger than the 
wavelength of the sine waves. The residual, b — B ~ B, is nor- 
mally referred to as the small-scale or fluctuating field, but in the 
present case such a characterization might be misleading, because 
such a field is quite regular and not actually fluctuating in the real 
sense of the word. Note in particular that the forcing function has a 
finite average, i.e. 



= Po(sinfcoz,0,0), 



(10) 



so the residual is / = -Fo(0, sin kox, sin koy). It turns out that also 
U and B point mainly in the x direction. Throughout this paper we 
denote the residuals by lower case characters. 



3 SIMULATION RESULTS 

Dynamo action is possible once the value of Rm exceeds a cer- 
tain critic al value of around 3; s ee Fig .fT] A similar curve was first 
shown by lOallowav & Proctoil d 19921) for the case of a prescribed 
flow U — Uo- For smaller values of Pm the growth rate is negative 
while for larger values it levels off at a value comparable to itofco. 



|http : / /pencil- code ■ googlecode . corri| 
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Figure 1. Dependence of the dynamo growth rate A on . Note that the 
critical value of i?m for dynamo action is around 3. For larger values of /?ni 
the growth rate levels off at a value around u^kQ. 



One may expect the Archontis flow to be a small-scale dy- 
namo, which means that the scale of the field would not exceed the 
scale of the flow, Lo- In order to check whether this flow can also 
generate fields on a scale larger than that of the flow we consider 
now also cases with L = 2Lo and 4I/o- In Fig.|2]we compare vi- 
sualizations of Bx and Bz for L/Lo — 1, 2, and 4. In the case 
L = Lo the magnetic field has a scale that is equal to that of the 
flow, but in the other cases the field breaks up into smaller scale 
contributions with a modulation in the y direction on the scale of 
the domain. In the latter case, the field on the s cale of the domain i s 
reminiscent to that found in helical turbulence ( lBrandenburd200 ih . 
but it is less dominant and less persistent than for L — Lq. This is 
mainly explained by a strong reduction of net cross helicity when 
the field breaks up into smaller-scale contributions. For these rea- 
sons we focus in the remainder of this paper on the case L = Lo, 
which is perhaps the simplest case known to produce net cross he- 
licity. 

Those dynamos produce large-scale fields, but they are not as 
prominent and persistent as in the case of large-scale dynamos that 
are driven by kinetic helicity. This is mainly because in the simu- 
lations with larger domains the cross helicity is strongly reduced 
once the magnetic field breaks up into smaller-scale fields. 

In Fig. [3] we show the evolution of the mean magnetic field, 
mean velocity and the small-scale cross helicity, he — {u ■ b), for 
a run with 7?m ~ 16 in a logarithmic scale; see Panel 1 and also 
the evolution of the magnetic energy compared to that of an dy- 
namo on a linear scale in panel 2. Time is normalized with respect 
to the microscopic visco-resistive time scale, (^fcg)"^. Given that 
the initial magnetic field is spatially random, it is first smoothened 
by resistive effects, leading to a short period where the magnetic 
energy decreases. Exponential growth occurs after about half a 
visco-resistive time, and then turns into a slow saturation phase af- 
ter about two visco-resistive times, which is best seen on a linear 
scale (lower panel of Fig. [3}. However, the late saturation behavior 
deviates from that of the dynamo, where the late evolution of 
the mean field is well described by a switch-on curve of the form 

S'' ~ l-exp(-At/r^), (11) 

where At = t — ts is the time after the end of the exponential 
growth phase alt = Ib and — {2riki)^^ is the large-scale resis- 




O 



Figure 2. Visuahzation of Bx and on the periphery of the domain in 
models with L/Lq = 1 (upper row), 2 (middle row), and 4 (lower row) for 
Rui = 13. In the case L = Lg the magnetic field has a scale that is equal 
to that of the flow, but in the other cases the field breaks up into smaller 
scale contributions with a modulation in the y direction on the scale of the 
domain. The coordinate directions are indicated in the lower left panel and 
the origin is indicated by O. 

tive time based on the wavenumber ki, which would be equal to fco 
in the present case. 

In Fig. m we demonstrate that the saturation time is essentially 
independent of the value of R^. Here, time is expressed in dynami- 
cal units by normalizing it in terms of the turnover time (uoko)^^ . 
The amplitude of the mean field increases mildly with i?^. A suit- 
able non-dimensional representation of the mean field is the quan- 
tity fiko/ hoPoUq. This number turns out to be of order unity 
and only weakly dependent on the value of i?^ for values between 
5 and 20. 

In the three cases displayed in Fig.|4]we have verified that the 
choice of averaging is unimportant. In other words, the results for 
yz and xz averages agree with those for the xy averages shown in 
Fig. |4] within 0.1-0.5 per cent. 



4 TURBULENT MAGNETIC DIFFUSIVITY 

In a number of circumstances it has been possible to character- 
ize the production of mean magnetic field in terms of a effect and 
turbulent magnetic diffusivity. Here, "turbulent" refers to the com- 
monly used name for transport coefficients describing the evolution 
of mean fields rather than a distinction between turbulent versus 
laminar flow properties. Both a effect and turbulent magnetic dif- 
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Figure 3. Saturation behavior for a run with = 16. The 

dotted line shows that the simple-minded helicity constraint formula 
does not describe the saturation of coirectly. The labels U, B, 
and he denote (jy^ko / Fq)'^/'^ , (B^ko/ ^ioPoFq)^''^ , and {u ■ b) 
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Figure 4. Comparison of the saturation behavior of for three different val- 
ues of Rf_i . 

fusivity have been determined also for other laminar flows such as 
the Roberts flow jBrandenburg. Radler & Schrinneill2008h . How- 
ever, such a description may not be applicable in the present case 
because of the possible presence of the additional Yoshizawa effect. 
Ignoring this complication for a moment, we can determine the otij 
and riij tensors in the relation 

{u X b)i = aijBj — rjijjj (12) 

using the test-field method dSchrinner etai]|2005l l2007t) . In this 
approach one solves an additional set of three-dimensional par- 



tial differential equations for vector fields 6'"', where the labels 
p = 1, 2 and q — 1,2 correspond to different pre-determined one- 
dimensional test fields B'"'. This leads to four vector equations for 
u X b^i that allow us to determine all components of Uij and rjij 
as functions of z and t. Owing to homogeneity and stationarity, it 
makes sense to present their average s over z and t . The test-field 
method has been criticized by Catta neo & Hughe i l|2009) on the 
grounds that the small-scale dynamo action would affect the re- 
sults. However, for magnetic Reynolds numbers of up to about 100 
the results of the test-field method have been proven to be consis- 
tent with results from direct simulations jMitra et al]|2009l) . 

The evolution equations for b^'' are derived by subtract- 
ing the mean-field evolution equation from the evolution equa- 
tion for B. These equations are distinct from the original induc- 
tion equation in that the curl of the resulting mean electromo- 
tive force is subtracted. This method has been successfully ap- 
plied to the kinematic case of weak magnetic fields in the pres- 
ence o f homogeneous turbulence either withou t shear jSur et al] 
|2008'; 'Brandenburg, Radler & Schrinner' '2008^ or with shear 
(Brandenburg 2005; Brandenburg et al. 2008a), as well as to the 
non-kinematic c ase with equipartition-strength dynamo-generated 
magn etic fields jBrandenburg et alj2008bljTilgner & Brandenburg! 
l2008l) . 

Using this method, it turns out that all components of aij van- 
ish within error bars, and that rjij has only diagonal components. 
However, as shown in Fig. |5] the 7722 component can be negative 
within a limited range of wavenumbers. (The fact that rjn 7^ 7722 
is not a priori surprising, because both U and B have only compo- 
nents in the x direction.) One of the two growth rates, 

= + '7ii)fco, A2 = -(r? + r?22)fco (13) 

is therefore positive. This suggests that there is the possibility 
of driving a dynamo by a negative t u rbulen t resis tivity effect 
dZheligovskv. Podvigina & Prised I2OOII : lUrpiij I2OO2I) . In such a 
case it is important to determine the wavenumber where the growth 
rate is largest. In our case, this happens for fc ~ fci (see lower of 
Fig.H). 

In the following we discard the possibility of dynamo activ- 
ity driven through a negative turbulent resistivity effect, because 
the test-field method ignores the presence of the Yoshizawa effect. 
Thus, we argue that equation l ll2t is an inadequate ansatz that re- 
sults in an apparent negative turbulent resistivity component. In the 
absence of a proper method for determining rjij we consider now 
a phenomenological description of the Yoshizawa effect using an 
isotropic turbulent resistivity, rjt- 



5 PHENOMENOLOGY 

The slow saturation process found here is reminiscent of the slow 
saturation process found for the dynamo, where net magnetic 
helicity is being produced on a resistive time scale. In the present 
case the magnetic helicity is essentially zero, but net cross helic- 
ity is being produced. Owing to the conservation of cross helicity, 
there is the possibility here too that full saturation requires a visco- 
resistive time scale, = (/ikf)^^, where /i — u + ri and fcf is the 
wavenumber corresponding to the typical scale of u and b. In our 
case, these fields depend essentially only on the x and y directions, 
so fcf = 2fcQ. The form of this relation is not known, although it 
is already clear that it is not the same as in the case of the nonlin- 
ear a effect. Most importantly, the saturation time does not seem 
to depend sensitively on the value of (Fig.|4j. Moreover, owing 
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Figure 5. Dependence of the normalized diagonal components of the tur- 
bulent resistivity tensor for iJni = 6 (upper panel) together with the coiTe- 
sponding growth rates (lower panel). 



to the presence of a forcing term in the momentum equation, the 
cross helicity is not necessarily conserved in the limit ^ Q, but 
it may change. Indeed, under the assumption of incompressibility, 
the evolution of the cross helicity per unit volume, (C7 ■ B) , is given 
by 



At 



(U-B)^ {F ■ B) - ii{W ■ J). 



(14) 



Here, angular brackets denote volume averages and W — 'V x U 
is the vorticity. Note the presence of the forcing term that can lead 
to the production of net cross helicity if the field has a component 
that is aligned with the forcing. 

Next, we restrict ourselves to horizontal averages, denoted by 
an overbar, and consider first their evolution equations. 



d_ 

dt 



A = UxB + £-r)J, 



^UxW + JxB + F + T-iyQ, 

at 



(15) 



(16) 



where £ = u x 6 is the mean electromotive force due to the 
correlation of small-scale velocity and magnetic field correlations, 
^ = uxw + jxb is the mean force due to advection and 
Lorentz force of small scale contributions, and Q = ^ x W is 
the curl of the vorticity. As discussed above, lower case characters 
denote the residual or "fluctuating" components, so for example 
w — W — W is the residual vorticity. 

We note that the U x W and J x B terms will be of no sig- 
nificance, because for our one-dimensional ^-dependent averages 
only the x and y components of A and U will be important for the 
evolution of the dynamo. We assume that £ has only contributions 
from thelVoshizawa ( 199G) effect and from turbulent resistivity and 
that J- has only a contribution from turbulent viscosity, i.e. 



£ = rW - r)tJ, 



(17) 



T=-vxQ. (18) 

A simplified derivation of the lYoshizawal ( fl990l) effect is given in 
Appendix [a] which shows that 



T — TU ■ b. 



(19) 



We have chosen here the symbol T instead of Yoshizawa's original 
symbol 7, because 7 is frequently used to describe the turbulent 
pumping velocity. Furthermore, T looks similar to 7 and it also 
reminds of the letter Y in Yoshizawa's name. 



In addition, there is also turbulent viscosity i/t 



_L5J 



and tu rbulent resistivity rjt = ^tu^ ( iKitchatinov. Riidiger & Pit^iii 
Il994h. although numerical simulati ons suggest Ut ~ 77* 
l lYousef. Brandenburg & Rtidigeill2()03h . Here, t is a typical time 
scale that may be estimated in terms of the turnover time, r — 

(«rnisfco)"\ where ttrnis = {u^)^^^. 

Inserting equations l ll7t and l llSt into equations il5\ and il6i . 
we derive the following evolution equation for the cross helicities 
of the mean and fluctuating fields: 



dt 
d_ 

dt 



{U-B) = {F-B) + r{W^) - fn{W ■ J), 



(20) 



{u-b) = (f-b)- r{W^) + ■ J) - ti{w ■ 3), (21) 



where = i/t+ Vt is the sum of turbulent viscosity and resistivity 
and ~ Ht + jJ. is the total (turbulent and microscopic) value. 
One can easily verify that the sum of equations l |20t and | |2U gives 
equation d. 

In the following we shall use equation l l21t to describe the 
evolution of T fully in terms of mean field quantities. This ap- 
proach was recently perused bv lKandusl(l2007i) for the more com- 
plete case where kinetic and magnetic helicities are also present. In 
equation i2H the term {u ■ b) is directly related to the mean field 
quantity T, and so is {w ■ j) — kf{u ■ b). An exception is the cor- 
relation of the forcing term with b, i.e. the term {f ■ b). However, 
it turns out that for the Archontis flow considered here, each of the 
three terms, (FiBi) for i = 1, 2, and 3 contribute equal amounts. 



so (F ■ B) = i(F ■ B) and (/ • b) 
express 

{f-b)^2{F-B) 



;(F ■ B),so that we can 



(22) 



purely in terms of mean field quantities. The validity of these re- 
lations can be seen in Fig. [6] where we plot the aforementioned 
correlations for a run with Rui = 32. 

With these preparations we can write down an evolution equa- 
tion for T, 



dT 



2t{F ■ B) ~ tTW^ + ntT{J ■ W) 



R, 



(23) 



At ■ r 

where we have defined a modified visco-resistive Reynolds number 



Note that it is related to R^^ via 

-Rm = {kf/kof{u^n^s/uo)Rfj.. 



(24) 



(25) 



Analogous to the magnetic case we can write this equation as a 
quenching formula by keeping the time derivative as an implicit 
term. 



r ^R, 



2t\F ■ B) + rVt(iy ■ J) - rAT/At 



(26) 
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Figure 6. Plot of the force-magnetic field correlations for a run with 
32. 



Figure 7. Saturation behavior of large-scale and small-scale cross helicities 
for two different values of iJm- 



These equations show that the generation of large-scale magnetic 
field by the T term produces (U ■ B) of the same sign as that of 
T (= TU ■ b). This is also seen in the simulations, where {U ■ B) 
and u • b do indeed have identical signs for the same Rui, but could 
individually, depending on initial conditions, have different signs; 
see Fig. |7] for two cases with different i?ni values. However, this 
sign property is a major difference to the case of the dynamo 
where {A ■ B) and (a ■ b) have opposite signs. The reason for this 
lies in the absence of a T term that is independent of u ■ b, i.e. there 
is only the term T = rit ■ t>. By contrast, the a effect has also a 
contribution from Idnetic helicity that is independent of magnetic 
helicity, i.e. a = ^rj ■ b — |rtj ■ u. 

It is instructive to inspect this difference by comparing equa- 
tion ( 1261 ) with the anal ogous equation for a quenching. Written in 
implicit form (see, e.g.. lBrandenburg|2008h , and ignoring magnetic 
helicity fluxes, this equation takes the form 



ao + R,^{rit{J ■ B)/B^^ 



Tcla/dt) 



(27) 



l + Rn,{B'')/Bl, 

^^^'^^ is the equipartition field strength and ao is 
the kinematic a effect, i.e. the term proportional to lj ■ it, which is 



where Bcq = {pu ) 



the crucial term that has no correspondence with equation l l26b . An- 
other difference is the presence of the forcing term in equation l l26b . 
Apart from that the two equations are quite analogous, i.e. is re- 
placed by R^, f^t is replaced by r/t, T is replaced a, t'^ is replaced 
by B^q , and W is replaced by B. 



6 KINEMATIC GROWTH PHASE 

The equations discussed above were originally motivated by try- 
ing to understand the late nonlinear stage of the dynamo. However, 
as we see from equation i26l . T itself has terms proportional to 
the mean field, suggesting that T should increase with the mean 
magnetic field. This is indeed the case during the kinematic stage; 
see the dotted line in the upper panel of Fig. |3] This suggests that 
the T term might also be responsible for the kinematic exponential 
growth of the dynamo. In order to identify the relative importance 
of this mechanism compared with the negative magnetic diffusiv- 
ity effect discussed at the end of Sect. [3] we investigate a simple 
model based on the induction and momentum equations along with 



the evolution equation for the small-scale cross helicity. The z- 
dependent averaging procedure for the mean magnetic and velocity 
fields then implies. 



dB 
'dt 



V X (TW -77tJ), 



(28) 



(29) 



where i^t = Vt + v and r]T ~ i]t + rj. We write Equations l |28l l 
and l |29t along with equation l l23t in the form, 

U^Fo- UTklU, (30) 

B ^TklU --qTklB, (31) 

t = 2tFoB + TklUiiitB - TU) - 7?~V"'T, (32) 

where the dots denote a time derivative, and double z derivatives 
have been replaced by a multiplication with — fcg. During the early 
kinematic phase the mean velocity is approximately constant. Sim- 
ulation results for a run with _Rm = 32 then yield U = U /uq ~ 
0.6. Applying therefore equation l |30l l to the steady state gives 
vi = Fo/koU, i.e. ut = 1.7uo/fco- 

The early exponential growth of both B and T is governed by 
just the first terms on the r.h.s. of equations (jSTJ and l |32| l, i.e. 



d_ f B 
dt [ T 





2rFo 



k^oU 




(33) 



This assumes that the i]t term in equation ( |31| l is negligible. There- 
fore the expected maximal growth rate for the Yoshizawa effect is 



(34) 



At = ±y/2FoTklU. 

Here we may estimate r in terms of the turnover time, r = 
(urmsfco)^^. Our dimensiouless turnover time, ito/itrms, is then 
about 0.4, so the dimensionless growth rate is 



uofco 



(35) 



This amounts to about 0.7, which is in good agreement with the 
simulation data. This suggests that the Yoshizawa effect may in- 
deed be responsible for driving the dynamo in the kinematic stage. 

This simple model does not describe the nonlinear satura- 
tion process. So, if one wanted to model this, one would need to 
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assume some ad hoc quenching prescriptions for various quanti- 
ties such as Vt, rjt, and r. This is in stark contrast to the case 
of the dynamo where equation l l27t describes both the kine- 
matic growth and the slow saturation phase quite accurately in 
the case of periodic boundary co nditions jpield & Blackmaiill2002l : 
iBlackman & Brandenbur3l2002h . 



7 CONCLUSIONS 

We considered here the Archontis flow, which is a generalization 
of the ABC flow. Such a flow was thought to be a small-scale dy- 
namo capable of generating magnetic fields at most on the scale 
of the flow. However, this flow tends to produce net cross helic- 
ity, whic h can lead to a m ean-field dynamo effect proposed orig- 
inally by I Yoshizawal ( ll99Q) . Direct numerical simulations of such 
flows performed with bigger box size show the presence of mag- 
netic fields on scales larger than the scale of the box (Fig.O. This 
is reminiscent of large-scale dynamos driven by kinetic helicity, 
where the resulting field is however much more prominent or per- 
sistent. 

The strongest cross-helicity production is found when the 
scale of the domain coincides with that of the flow. In that case 
dynamo action is possible once Rm exceeds a certain critical value 
which in our units turns out to be Rm — 3. The present work has 
shown that the kinematic phase of the Archontis dynamo can be 
modelled in terms of the Yoshizawa effect. The sign of the cross 
helicity depends on initial conditions, so either sign is possible for 
one and the same flow field. Simple phenomenological considera- 
tions support the idea that the Yoshizawa effect can be expressed 
in terms of the mean field alone; see equation J26b . This expression 
looks similar to the dynamical quenching formula for the a effect 
under the constraint of magnetic helicity conservation. However, 
this expression does not actually describe quenching, but growth. 
So, contrary to our initial expectation, this mechanism is not con- 
stant in time and so it does not correspond to a battery with linear 
growth, as was assumed by Brandenburg & Urpin ( 1998). Instead, 
it leads to exponential growth. At the end of the exponential growth 
phase the dynamo shows a characteristic saturation behavior that is 
reminiscent of a effect dynamos that are controlled by resistive 
magnetic helicity evolution. In the present case, the conservation 
of cross helicity was initially thought to be responsible for this pro- 
longed saturation behavior, but it turns out that the presence of a 
forcing term in the momentum equation can lead to a production of 
net cross helicity even in the ideal limit. 

It has long been speculated that the Yoshizawa effect could 
be relevant in accretion disc s and galaxies where differential rota- 
tion is strong ( lYokoilll996h . However, it turns out that, unlike a 
effect dynamos that normally have a given kinematic value of a, 
the T term cannot be calculated a priori, but it itself depends on the 
mean field. The end result is again reminiscent of the a effect in 
that both TW as well as aB are linear in B during the kinematic 
growth phase. However, the results of the test-field method show 
clearly that there is no a effect in that case. Indeed, the mean elec- 
tromotive force has no component along the mean magnetic field, 
confirming that the re is no a effect. There is also no shea r-current 
(or Wx7) effect feoga chevskii & Kleeorinl2003Ll2004h . because 
the off-diagonal components of rjij were found to be zero within 
error bars (Sect. I?^. This supports the idea that the growth of the 
magnetic field is here indeed the result of the Yoshizawa effect. Al- 
though the 7722 component is found to be negative when ignoring 



the Yoshizawa effect, it is argued that this result is an artifact of 
using an inadequate ansatz for the mean electromotive force. 

Obviously, the flow considered here is relatively simple and 
hardly of direct astrophysical relevance. However, it has been sug- 
gested that dynamos with field-aligned flows might be particu- 
larly efficient in generating magnetic fields in the solar tachocline 
(Galloway 2008). If those ideas can be substantiated, it would be 
interesting to see whether the phenomenological description devel- 
oped in the present paper carries over also to other cases such as 
this tachocline model. 

Another possible avenue for future research would be the 
study of fully turbulent dynamos in the presence of cross helicity. 
An example of this was shown in Fig. |2] where the flow was driven 
by the Archontis forcing function, but on a scale that is smaller than 
that of the computational domain. Those dynamos produce large- 
scale fields, but they are not as prominent and persistent as in the 
case of large-scale dynamos that are driven by kinetic helicity. This 
is mainly because in the simulations with larger domains the cross 
helicity is strongly reduced once the magnetic field breaks up into 
smaller-scale fields. 

As mentioned in the introduction, the solar wind is one of the 
few examples where the turbulence is believed to have net cross 
helicity, but with opposite signs in the two hemispheres. Although 
the solar wind is not normally thought to harbor dynamos, there 
is the problem of an unexplained contribution to energy deposition 
away from the source. It would therefore be worthwhile exploring 
the role of the Yoshizawa effect in the conversion of energy in solar 
wind turbulence. 
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where we have included only the antisymmetric contribution to 
VL/" that leads to terms with W, i.e. Uij — —^>^ijkWk+ the 
symmetric part, where a comma denotes a partial derivative. Next, 
we calculate d£ /dt and include only terms proportional to it • 6 by 
assuming Uibj = ^SijU- b-\- terms proportional to u x b, but those 
would later not contribute to the component of £ that is parallel to 

W . In this way we obtain from u x b and ii x b each the term 
|w ■ b, so 

d£/dt — I It • 6 + ... — triple correlations. (A4) 

In the spirit of the minimal r approximation we approximate the 
triple correlations by a quadratic correlation in the form of a damp- 
ing term, i.e. we assume that the triple correlations are equal to 
£/t. Finally, assuming stationarity, we drop the time derivative and 
obtain £ = |tw ■ b. 

In an alternative derivation one can write U-'VU = U x W — 
|VC/^ and subsume the gradient term in a generalized pressure 
term. Splitting U x W into mean and fluctuating part yields then 
directly a term u x W without the 1/2 factor The final result is 
then 

£ = TiTb, (A5) 
which is also the expression used here. 



APPENDIX A: CROSS-HELICITY EFFECT 

We present here a simplified derivation of equation i26\ using the 
minimal r approximation. We use the linearized evolution equa- 
tions for the fluctuations b and u and calculate 

d£/dt = u xb + uxb. (Al) 

In order to highlight the essence of the lYoshizawl (Il990h term we 
isolate from the very beginning the terms that are proportional to 
the mean vorticity. Thus, we consider in the evolution equations of 
b and u only those terms that contribute to terms proportional to 
W and write 

b^+b-VU + ... = -\bxW+..., (A2) 
u = -u-VTJ + ... = +luxW + (A3) 



© 2009 RAS, MNRAS 000,[T]{8] 



